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Collisions of deformed uranium nuclei provide a unique opportunity to study the spatial de- 
pendence of charmonium in-medium effects. By selecting the orientations of the colliding nuclei, 
different path lengths through the nuclear medium could be selected within the same experimental 
environment. In addition, higher energy densities can be achieved in U+U collisions relative to 
Au+Au collisions. In this paper, we investigate the prospects for charmonium studies with U+U 
collisions. We discuss the effects of shadowing and nuclear absorption on the J/ij) yield. We intro- 
duce a new observable which could help distinguish between different types of J/ip interactions in 
hot and dense matter. 

I. INTRODUCTION 

Collisions of deformed nuclei, such as 238 U, are an interesting alternative to those of more spherical 197 Au nuclei 
at RHIC because changes in the nuclear orientation allow wider variations in energy density within the same system 
PHI]. Previously, studies of U+U collisions were discussed in the context of elliptic flow. The quantity V2/e, where V2 
is the elliptic flow parameter and e is the initial spatial eccentricity, provides valuable information about the matter 
created in heavy-ion collisions [5]- In a dilute system, w 2 /e scales like 1/ S±dN c h/dy, where dN^/dy is the charged 
particle density and S± is the transverse area of the overlap zone. In central Au+Au collisions at y/s NN — 200 GeV 
i>2/e reaches the limit predicted by ideal hydrodynamics for strongly-interacting matter [6] which can be interpreted 
as a signature of quark-gluon plasma (QGP) formation. It was argued that U+U interactions could provide higher 
densities and a larger spatial eccentricity, testing whether vijt saturates, as predicted by ideal hydrodynamics, or 
increases further with dN^/dy [TJ|4l[7]. 

Charmonium production is another key observable for studying the properties of the hot and dense matter created 
in relativistic heavy-ion collisions. Some time ago J/ip suppression was proposed as a signal of QGP formation [8]. 
The suppression is expected to arise from color screening of the binding potential in a QGP, similar to Debye screening 
in a classical electromagnetic plasma. The magnitude of the suppression depends on the charmonium binding energy 
and the energy density of the medium, related to the temperature. Therefore, studies of J/ip production, particularly 
at low transverse momenta, reveal the thermodynamic properties of the medium. 

The suppression of J/ip production has been studied in detail at the CERN SPS and at RHIC. The NA50 and 
NA60 experiments at the CERN SPS, ^/s NN = 17.3 GeV, observed strong J ftp suppression as a function of collision 
centrality ffiUlj. Results from the PHENIX collaboration at RHIC show that the J/ip suppression in midrapidity 
Au+Au collisions at y / s NN = 200 GeV is similar to that observed at SPS energies [T2] even though the energy 
densities and temperatures reached at RHIC are much higher than those at the SPS. Moreover, the forward rapidity 
J/ip suppression at RHIC is stronger than at midrapidity. Such a pattern suggests that either additional processes 
compensate for the effects of color screening or the suppression has different origins at the two energies. 

There are several processes which may affect charmonium production in nucleus-nucleus collisions. For comprehen- 
sive reviews of different aspects of quarkonium in-medium interactions, see Refs. |13l 114] . Here we briefly summarize 
a few relevant features, some of which we address in more detail later. 

Suppression by cold nuclear matter (CNM) effects such as parton shadowing and absorption by nuclear matter are 
very important, although not yet established with satisfactory precision [15] . For example, parameterizations of gluon 
shadowing differ significantly between models [16] . At RHIC energies, CNM effects are studied in d+Au collisions. 
Changing the shadowing parametrization can significantly affect the magnitude of the absorption cross section required 
to match the measured nuclear modification factor i?dAu- I n addition, feed-down from radiative or hadronic decays 
account for 40-50% of the observed 15 quarkonium states (J/ip and T). The higher mass quarkonium states have 
very different radii and formation times and should thus also have different absorption cross sections. While the 
experimentally measured ip' absorption cross section is significantly larger than the J/ip at fixed-target energies, the 
difference decreases with increasing ^/s NN |171 118] . The energy and A dependence of \c absorption is unknown |19) . 
The importance of absorption for higher bottomonium states is also not known. 
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If color screening dominates the J/ip in- medium interaction, the characteristic pattern of sequential charmonium 
"melting" should be observed. The binding energies of higher mass charmonium states, ip' and \c, are smaller than the 
J/ip. Therefore these states dissociate at lower energy densities compared to J/ip so that, at higher T their feed-down 
would no longer contribute to the measured J/ip yield. The remaining, direct, J/ip production would be suppressed 
at larger energy densities. This process would result in a step-like dependence of the J/ip survival probability as a 
function of energy density [20] , A similar sequential melting would apply to bottomonium. 

Suppression due to interactions with comoving partonic and hadronic matter has also been postulated. The A 
dependence of comover dissociation is similar to that of nuclear absorption [211 122) . The effects are generally assumed 
to be small although they strongly depend on the model [2TJ [23J [53] . 

J/ip production may be also enhanced by statistical coalescence in QGP. In this case, the J/ip yield, A^{f , is 
proportional to the square of the number of charm quarks in the system, N c , so that N^/Jf oc N 2 . The charm yield 
is expected to be proportional to the number of binary NN collisions, N hin , i.e. AT cg oc iV bin . The cc production 
rate is not expected to be modified in the QGP. If the density of charm quarks is high enough, then secondary J/ip 
production by coalescence of uncorrelated c and c quarks at hadronization can occur |25l I26j . 

Combinations of the above effects have previously been used to describe the J/ip data. Despite different assumptions 
about the initial state of the system, models including cold matter effects together with different combinations of color 
screening, comover absorption, and coalescence give quantitatively similar results at RHIC [27) . It is expected that 
U+U collisions could provide additional means for distinguishing between these effects. 

By selecting particular orientations of U+U collisions, different path lengths through nuclear matter, resulting in 
different effective nuclear absorption cross sections and available energy densities, could be studied. Interactions with 
the long axes of U nuclei aligned along the beam ( "Tip+Tip" or "TT" configurations) give the highest energy densities 
as well as the longest path, L, through the matter. In contrast, L is shortest in the configuration where the short axes 
of the U nuclei are aligned along the beam axis ("Side+Side" or "SS" configurations). Different orientations of U+U 
interactions could study nuclear absorption with reduced uncertainties due to shadowing, as we will show. In general, 
U+U interactions provide an additional check on models which describe the Au+Au data because some effects, such 
as shadowing, should be similar in U+U and Au+Au interactions. 

We investigate the feasibility of complementing Au+Au collisions with U+U studies of charmonium production and 
suppression to distinguish between scenarios of J/ip in-medium interactions. We also discuss CNM effects in U+U 
collisions. 



Here, we discuss two types of cold nuclear matter effects: nuclear absorption and shadowing. Shadowing, the 
modification of the parton densities in nuclei relative to free protons, is an initial-state effect. Nuclear absorption, the 
breakup of the charmonium state as it traverses nuclear matter, is a final-state process. The two effects are assumed 
to factorize. We neglect interactions with hadronic comovers. 

We employ the computational framework for studies of cold nuclear matter effects described in Ref. [28] . The 
nuclear parton densities, f^(x,Q 2 , fy, z), where A is the atomic mass number, ?t and z are the transverse and 
longitudinal location of the parton in position space, x is the parton momentum fraction and Q 2 is the interaction 
scale, factorize into the nucleon density in the nucleus, pA{fr,z), independent of the kinematics; the nucleon parton 
density, f[(x,Q 2 ), independent of A; and a shadowing ratio, S' l (A, x, Q 2 , rV, z) that parameterizes the modifications 
of the nucleon parton densities in the nucleus so that 



The effect of shadowing is stronger in central collisions, while at asymptotic distances shadowing disappears. Averaging 
over impact parameter, 6, the minimum bias result measured in nuclear deep-inelastic scattering is regained. We 
use the EPS09 shadowing parametrization [16] and assume that the impact parameter dependence of shadowing is 
proportional to the local nuclear density [29j and exhibits a relatively sharp transition region around b ~ 2Ra, see 
Fig.[l] If one assumes instead that the impact parameter dependence is proportional to the path length L, S' 1 exhibits 
a weaker transition region [30] . We discuss the consequences of different assumptions of this b dependence in the next 
section. At RHIC energies quarkonium production is dominated by gluon fusion. Therefore, shadowing generally 
refers to gluon shadowing in this context. We employed the Woods-Saxon density distribution for deformed nuclei [3] 
to describe pxj(rr,z), 



II. COMPUTATIONAL FRAMEWORK 



tf(x,Q 2 ,r T ,z) = PA (f T ,z)S i (A,x,Q 2 ,f T ,z)ff(x,Q 2 ) . 



(1) 



pu(r T ,z) = p (l+cxp[(r - R c )/a}) 1 



(2) 
(3) 
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where po is the density at the center of the nucleus, Rc is the total effective radius, Y20 and I 40 are Leg endre 
polynomials, dependent on cos 9, where 6 is the polar angle, fx — (?t cos#, tt smO) and r = a/I^tI 2 + z 2 . The 
deformation parameters are f3 2 — 0.28 and /?4 = 0.093 with R — 6.81 fm and a = 0.54 fm. In the SS orientation, 
cos0 = and Rc = Rs — 6.41 fm while in the TT orientation, cosd = 1 and Rc — Rt — 8.56 fm. 

To implement nuclear absorption on quarkonium production, the production cross section is weighted by the survival 
probability, S abs , 

S ahs {b-r T ,z') =cxpj-^ dz"p A {b-r T ,z")a ahs {z" -z')} (4) 

where z' is the longitudinal production point and z" is the point at which the state is absorbed. The nucleon absorption 
cross section, <7 a bs, typically depends on the spatial location at which the state is produced and how far it travels 
through the medium. 



III. COLD NUCLEAR MATTER EFFECTS IN U+U COLLISIONS 



In this section, we investigate whether it is possible to separate the CNM effects of shadowing and absorption in U+U 
collisions. Nuclear absorption depends on the nuclear path length L. Therefore, it is expected to be considerably 
different in the TT and SS orientations of U+U collisions. Because L also depends on the impact parameter, b, 
the effective magnitude of nuclear absorption also changes with b. The dependence of shadowing on the TT and SS 
orientations is expected be less significant. Moreover, if the difference between collisions in the TT and SS orientations 
does not depend (or depends only weakly) on the strength of the shadowing, then the ratio of quarkonium production 
in TT and SS orientations would be effectively independent of the shadowing parametrization. This would greatly 
simplify the extraction of absorption effects because the systematic error due to shadowing would be significantly 
reduced. 

We first investigate shadowing and nuclear absorption separately and then discuss the combined intensity of these 
CNM effects in U+U collisions. In order to estimate the strength of each effect, we compare results in the TT and SS 
orientations where the difference between the two cases should be most apparent. We assume that it will be possible 
to select these configurations with an efficiency appropriate for charmonium studies, as discussed later. 

We begin with an investigation of the charmonium survival probability in U+U collisions for shadowing alone, S'shad- 
Figure [l] shows S s had as a function of b in TT and SS configurations of U+U collisions at y = and 1.75, calculated 
in the color evaporation model of quarkonium production 27] . (We do not show results at backward rapidity since 
the results are symmetric around y — in AA collisions.) Note that S s h a d < 1 represents shadowing, S s h a d > 1 
antishadowing, and S'shad = 1 no shadowing. The convolution of two shadowing parameterizations results in stronger 
shadowing at midrapidity than at forward rapidity for central EPS09 shadowing, as we show in more detail later. The 
CEM results exhibit a consistent pattern between shadowing effects at leading and next-to-leading order [15] , This 
behavior is typical of most shadowing parameterizations based on collinear factorization |31j . 

Calculations of color singlet J/?p interactions in the dipole approximation suggest that collinear factorization may 
be inapplicable due to the coherence of the interaction (35J|33]. In this case, higher- twist effects enhanced by powers 
of A 1 / 3 would dominate pA interactions and enhanced J/tp suppression should set in at large rapidity. However, such 
enhanced effects are seen in fixed-target energies as well, outside the range of validity of gluon saturation models 
[341 135) . and appear to scale with projectile momentum fractions, x\, rather than with the target fraction x 2 . 

While the most recent PHENIX data [36] seem to exhibit a stronger than linear impact parameter dependence 
of shadowing [37], the min bias results for i?dAu are in good agreement with calculations based on Ref. [15] . as are 
the results in the most central impact parameter bin. Indeed, these new data suggest that while the shadowing may 
decrease more strongly than the calculations of Refs. [33] [3U] predict, the strength increases rather slowly for low 
impact parameters. 

To illustrate the effect of a larger than linear impact parameter dependence of shadowing, we compare the results for 
shadowing parameterizations that depend on the local nuclear matter density and on the parton path length through 
nuclear matter 30 . 

S p (A,x,Q 2 ,f T ,z) = I + N p (S*(A,x,Q 2 ) - l) ^^^ , (5) 



S l TA (A,x,Q 2 ,f T ,z) = 1 + N Ta (S*(A,x,"- 



f J dzp A (r T ,z) \ n 
\ Jdzp A (0,z) J ' 



= l + N TA (S\A,x,Q 2 )-l)(^^y , i+) 
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Figure 1: (Color online) The survival probability due to shadowing, SWad in TT (lines) and SS (symbols) orientations of U+U 
collisions for (a),(c) y = and (b),(d) y = 1.75 assuming the impact parameter dependence of shadowing is proportional to pa 
(top) and Ta (bottom) for n = 1, 2 and 5. 



where Ta is the nuclear profile function and N p , Nt a are chosen to normalize the integral of S l (A, x, Q 2 , fr, z) over 
the nuclear volume, weighted by the nuclear density distribution, to S l (A, x, Q 2 ) in minbias collisions. The values of 
the normalizations are dependent on the power of n chosen as well as the parameterization of the impact parameter 
dependence. Results are shown in Fig. [T] for n = 1, 2 and 5. 

When the shadowing is parameterized according to the nuclear density, as in Eq. ^ , the difference between the TT 
and SS configurations is rather distinct over a broad range of impact parameters. The shape of S'shad(fr) resembles the 
Woods-Saxon density distribution. The normalization, N p , changes slowly with n. The difference between the two 
spatial orientations of U+U collisions, TT and SS, is very small for most central collisions, b < 6 fm. As n increases, 
N p increases, as does the difference between the TT and SS orientations. Nonetheless, the results retain the general 
Woods-Saxon shape while impact parameter dependence steepens over the transition from shadowing to free nucleon 
behavior. This transition is centered around b ~ 2Rq for n = 1, the impact parameter where the effect of shadowing 
is half the 6 = value. Note, that Rc = Rt, the length of the long axis in TT configurations, and Rc — Rs in SS 
configurations. 

When the longitudinal direction is integrated over, as in Eq. (J6j|, the distinctive Woods-Saxon shape seen in the 
upper half of Fig. [T] is somewhat washed out and Sghad increases rather smoothly with impact parameter. A separation 
between the TT and SS configurations is already apparent for small values of b. In this case, Nt a depends more 
strongly on n so that the difference between n = 1 and n — 2 is larger here than for Eq. d5D. 

Increasing n reduces the effective radius for shadowing in both formulations in Eqs. |5|) and j6j. We have only 
shown results for n < 5 but have checked higher values of n and see that the trend continues. The separation between 
the TT and SS configurations for Eq. ^ increases to 10% at low b for n — 15 while the impact parameter dependence 
of iSghad more closely resembles a step function. The growth of A^ p with n remains slow. On the other hand, with 
the dependence of Eq. ([6]), the difference between the orientations remains small, even at higher n, while Nt a is a 
stronger function of n. Large values of n are still consistent with the PHENIX data at forward rapidity [55]. 
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Figure 2: (Color online) The ratio of charmonium survival probabilities, Sshad, in TT and SS orientations assuming the impact 
parameter dependence of shadowing is proportional to pa (left) and Ta (right) for n = 1, 2 and 5. 
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Figure 3: (Color online) (a): The J/ip nuclear absorption survival probability, S a bs, in TT and SS orientations of U+U collisions 
for (Tabs — 3 and 5 mb. (b): The ratio of results in TT and SS configurations for cr a bs = 3 (solid) and 5 (dashed) mb. 



We expect that the impact parameter dependence can be characterized by other means and only show the possible 
range of effects here to study the sensitivity of CNM effects in U+U collisions on the impact parameter dependence 
of shadowing. A full analysis of the b dependence, under study in Ref. [38], is beyond the scope of this work. 

The ratio of S^ad m the two configurations is shown in Fig. [2] for the same values of y and n. Away from the 
transition region, the ratio is relatively constant and close to unity with a maximum 5% difference employing Eq. ^ 
with n = 5. The ratio of results with Eq. Q is much smaller at b ~ 0. Thus, for moderate shadowing, the ratio in 
central collisions is relatively insensitive to 5 s had- The largest change in the ratio with n = 1 is at ~ (Rt +Rs)/2. As 
n increases, the dip in the ratio TT/SS deepens and shifts to lower impact parameters. The dip is deeper for y = 
where the midrapidity shadowing effect is stronger. The forward rapidity effect is weaker but tends to shift to still 
lower impact parameters because the assumption of factorization pairs strong shadowing in one nucleus with small 
shadowing or some antishadowing in the other. In both cases, since f^(x,Q 2 ,rT,z) is proportional to the nuclear 
density distribution, see Eq. ([T]), the shape of .Sshad(fr) directly reflects the spatial orientations of the U nuclei. 

We now discuss the impact parameter dependence of the survival probability for nuclear absorption. Figure [3] (a) 
shows S^bs as a function of b for U+U collisions in the TT and SS configurations assuming Cabs = 3 and 5 mb. 
These two values of dabs bracket the range obtained by the RHIC experiments. As previously discussed, the min bias 
result, i?dAu, reported recently by PHENIX |3B] is relatively well described by the EPS09 shadowing parameterization 
with Cabs = 4 mb. Before the most recent PHENIX data [33] were available, the divergence of Rcp{y) from the 
calculations of absorption with shadowing was quantified the extraction of a rapidity-dependent absorption cross 
section |35j . These results showed a strong increase of the effective absorption cross section at forward rapidity which 
could be attributed to a heretofore neglected effect such as initial-state energy loss. However, a reanalysis of these 
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data based on a more complete understanding of the impact parameter dependence of nuclear shadowing would affect 
the magnitude of the extracted absorption cross section. This will be addressed in Ref. [38] . 

While a similar trend as a function of b is observed relative to Fig. [T] the difference between TT and SS orientations 
at b = is larger for absorption. Thus nuclear absorption is the only relevant orientation-dependent CNM effect in 
central collisions. The larger effect is due to the exponential factor in S'abs, proportional to the path length, in addition 
to the overall dependence on the density in the calculation of the total J/ip yield. Thus the angular orientation of 
the colliding nuclei affects absorption more strongly than shadowing. The exponential dependence on density in S'abs 
also broadens the transition region around b ~ 2Rc- 

The ratio of survival probabilities for the two orientations is shown in Fig. ^b) for er a b s = 3 and 5 mb. The ratio 
shows a ~ 20% difference relative to the maximum 5% for shadowing alone. It also varies more slowly with impact 
parameter, as seen by comparing Figs. [2] and [3]jb) . 

Figure [4] shows the ratios of the total J/tp survival probability in the TT and SS configurations due to cold nuclear 
matter effects, Sj/^ — S a \, s S s ha.Ai & s a function of impact parameter. In central collisions, b < 5 fm, the magnitude 
of (S TT /S ss )j/^ is most sensitive to cr a b s while at larger impact parameters, b > 10 fm, (S TT /S ss )j/^, is dominated 
by shadowing. If the TT and SS orientations could be effectively selected in an experiment, these features could help 
differentiate cold nuclear matter effects from color screening effects. However, the best case scenario would be to first 
attempt to select the tip and side orientations of the uranium nuclei in d+U collisions. 
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Figure 4: (Color online) The ratio of the cold nuclear matter J/tp survival probability in TT relative to SS orientations of U+U 
collisions for <j a bs = 3 and 5 mb assuming different parameterizations of the impact parameter dependence of shadowing (see 
text for details). 

Even if the TT or SS configurations cannot be selected due to either technical difficulties or low efficiencies, U+U 
collisions averaged over all orientations could be used for quarkonium studies. There are no significant experimental 
difficulties and the nuclear path length L and charged track density are respectively 5-10% and 10-15% larger for 
the same centrality class in Au+Au collisions [3] . The larger L increases the effective nuclear absorption. Figure [5] 
shows the orientation-integrated 5 a bs as a function of impact parameter for Au+Au and U+U collisions with cr a b s = 
3 and 5 mb. The difference between S^ b v s and 5 , ^ Au is largest, ~ 5%, at b ~ 11 — 12 fm and can be attributed to the 
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difference in radii, 6.8f fm for uranium and 6.38 fm for gold. The magnitude of the difference depends on cr a bs- 

On the other hand, the relative difference in the shadowing effect is negligible. The expected S s had is shown in 
Fig. [6] for d+Au and d+U collisions (a) as well as Au+Au and U+U collisions (b). Shadowing effects in U+U and 
Au+Au collisions at the same ^/s NN are very similar at a given value of b. Thus the uncertainty due to shadowing 
will cancel in the ratio of the J/ip yields in U+U and Au+Au collisions. 

We now discuss the modification of T production in Au+Au and U+U collisions due to cold nuclear matter effects. 
The nuclear absorption of Y(IS') seems to be significantly weaker than for J/ip |28) . Moreover, antishadowing may be 
expected at midrapidity rather than shadowing, supported by the experimental results at RHIC. The cross section of 
T production in d+Au collisions is well described by the color evaporation model with antishadowing and no nuclear 
absorption [40]. The expected Sshad for T is shown in Fig. [7] for d+Au and d+U collisions (a) as well as for Au+Au 
and U+U collisions (b). The effect of shadowing is much smaller than for J ftp. Instead, antishadowing at midrapidity 
may increase the T rate by as much as 20-30%. Note that antishadowing effects in U+U and Au+Au collisions at 
the same t/s~^/ are very similar, as is the case for the J/ip. Furthermore, a different formulation of the quarkonium 
production model, such as in Ref. (41], would yield similar results even though the larger average scale used in those 
calculations further reduces the overall shadowing effect. 




b [fm] b [fm] 

Figure 5: (Color online) Nuclear absorption as a function of b for Au+Au and U+U collisions averaged over all orientations 
(a) and the ratio of U+U relative Au+Au collisions (b). 

Finally, we comment on the suggestion that double-color filtering and mutual boosting of the saturation scales in 
colliding nuclei make the transition of cold nuclear matter effects from pA to AA collisions nontrivial [33J . However 
these effects are rather small at RHIC energies (less than 10% for double-color filtering and 18% for mutual boosting 
of the saturation scales). In addition, they act in the opposite directions, therefore canceling to a large extent. 
Reference [33J proposed that an increase of J/ip pr broadening in AA collisions relative to pA would directly reflect 
the boosting effect. Such an increase is not seen at RHIC. In fact, (p\) at midrapidity in Au+Au collisions at ^s NN 
— 200 GeV is approximately independent of centrality [39) . Determining the nature and strength of such a boosting 
effect at RHIC requires high precision data. Moreover, the effects proposed in Ref. [33] will be very similar in Au+Au 
and U+U collisions and thus cancel if ratios of J/ip production in the two systems are studied. 

IV. ENERGY DENSITY AND TEMPERATURE IN U+U COLLISIONS 

We now investigate what energy densities and temperatures are accessible in U+U collisions. The charged track 
density is expected to increase in U+U interactions compared to Au+Au. Therefore, the U+U temperature should 
be higher in the deconfined phase. Employing U+U collisions could extend the range of temperatures accessible at 
RHIC and test the sequential melting hypothesis of charmonium and bottomonium states. Suppression of T(IS') 
beyond CNM effects would constitute direct proof of color screening since secondary T production by coalescence 
from a QGP is negligible at RHIC. Comover absorption is also insignificant. Recent STAR results show that the 
T(IS') is not strongly suppressed in central Au+Au collisions [32]. Below, we investigate whether the temperature of 
the matter produced in U+U collisions could be high enough to expect "melting" of the Y(1S) state in addition to 
J/ip suppression. 

We use the Bjorken formula [33J to estimate the energy density available in U+U collisions. The Bjorken energy 
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Figure 6: (Color online) The J/ip survival probability due to shadowing in d+A (a) and A + A (b) collisions. The solid lines 
(filled circles) represent results for Au (U) nuclei obtained with the central EPS09 shadowing parametrization while the dashed 
lines (open circles) show the results with the upper and lower limits of the EPS09 uncertainty. 
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Figure 7: (Color online) The T survival probability due to shadowing in d+A (a) and A + A (b) collisions. The solid lines 
(filled circles) represent results for Au (U) nuclei obtained with the central EPS09 shadowing parametrization while the dashed 
lines (open circles) show the results with the upper and lower limits of the EPS09 uncertainty. 



density in a one-dimensional longitudinally expanding system can be calculated using 

1 dE T 



tS± dy 



(7) 



where r is the formation time of the medium, S± is the transverse area overlap of the colliding nuclei, and dE? /dy 
is the transverse energy density. We estimate dEx/dy from PHENIX data [H], 



dEr = E T dN ch 

dy v ^ v N ch dr, 



(8) 



where Et/N^ = (dE? / drf) / (dN^/ drf) is the average transverse energy density per charged track and C v ^ y = 
1.25±0.05 is a scale factor converting dET/dij\ v= Q to dEx/dy\ y —Q. In Au+Au collisions at y/s NN = 200 GeV, Et/N^ 
is constant over a broad range of event centralities (0-60%). We assume that (Et / W c h) AuAu is a good approximation 
of {E T /N ch ) vlJ . 

The most important ingredient is the overlap area S±, calculated in the Glauber framework. It depends strongly 
on the definition used in calculations. For example, the values of S± published by PHENIX [44] and STAR [45] 
are a factor of five larger than results calculated for U+U and Au+Au collisions in Refs. [21 E]. When the relative 
differences between the energy densities of U+U and Au+Au collisions are considered, this discrepancy is unimportant 
as long as a consistent definition is used. However it is important when the absolute value of ee is calculated. In this 
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paper, we consider the case where S± is defined as the transverse area of the overlap zone weighted by the number 
of participants [2j [3] . This definition gives the effective "hot" transverse area of the overlap zone but neglects areas 
with low participant density, a rather small effect. 

Figure [8] shows an example of the participant density profile for U+U collisions at b = in the TT and SS 
configurations as well as averaged over configurations calculated in the Glauber framework of Ref. [3] . Much higher 
density is observed in TT compared to SS configurations as well as relative to orientation-averaged U+U collisions since 
the average is dominated by SS configurations [46^ . Note also the narrower transverse profile of the TT configurations 
in Fig. |H 



(a) U+U Tip+Tip (b) U+U Side+Side (c) U+U Averaged over configurations 
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Figure 8: (Color online) Participant density profiles for U+U collisions in the TT and SS configurations as well as orientation- 
averaged collisions at b — 0. The circles are shown to guide the eye and are common to all plots [3l I46| . 



We use the results of Ref. [3] to estimate the increase in energy density in U+U relative to Au+Au collisions. We 
first estimate cb in orientation-averaged U+U collisions. Figure [9] shows the ratio of the product e^r for orientation- 
averaged Au+Au and U+U collisions. We note that es is 15 — 20% larger in U+U collisions relative to Au+Au 
collisions. Moreover, U+U collisions in the TT configuration could provide an increase of up to 30% in the charged 
track density, l/S±dN c \ 1 /dr), relative to orientation-averaged collisions at the same value of b, as shown in Fig. 10 'a). 
Since 1/ S±dN c h/di] is proportional to eg, TT configurations can thus increase egT by 20-30% in central and semi- 
central (b < 8 fm) U+U collisions, as shown in Fig. 10 b) . 
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Figure 9: The ratio of Bjorken energy densities in U+U and Au+Au collisions as a function of centrality. The shaded boxes 
are the systematic uncertainties on dN^ /dy and dN^ Au / dy added in quadrature. The statistical uncertainties are negligible 



In our case, the relevant formation time, r, is the time required to develop the qaurkonia wavefunction. The 
lower limit on the charmonium formation time is the time required for color neutralization of the cc pair, ~ 
(2m c AqcD) 1 ^ 2 — 0.25 fm for low px J/tp [17. . The formation time of the final-state J/V, however, may be somewhat 
longer. In Ref. |14j , the formation time for low pt ground state quarkonia is estimated to be the inverse of the binding 



energy, E bi nd, t 



1/E } 



bind • 



With R 



bind 



0.64 GeV and E^ ind 



1.1 GeV 



the formation times for the J/ij) 



10 



50 

M 45 

X 40 
:p 35 



03 



0, 



(a) 



° 30f—m 

25 
20 
15 
10 
5 



I TT 

► Average 
E - ▼ SS 



'I! 



Hi 



o 



10 



12 



14 
b [fm] 



3 « 1.4. 

« 1 2 
H 

h m 1 



I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 

■ — — ■ 




0.6 
0.4 
0.2 

0, 



►•••••••••••••••••••••**■■■ 



(b) 



TT/ Average 
SS/ Average 







10 



12 14 
b [fm] 



Figure 10: (Color online) (a): The average charged track density per unit of pseudorapidity 77 in U+U collisions in TT and SS 
configurations as well as orientation-averaged collisions as a function of impact parameter, (b): The ratios of energy densities in 
TT and SS orientations relative to orientation-averaged collisions as a function of b. In both cases, the statistical uncertainties 
are less than the size of the symbols. 



and T(1S) are tj/^ ~ 0.3 fm and tt ~ 0.18 fm respectively. We used two values of r in our studies: r = 0.25 fm 
to evaluate the maximum temperature relevant to J/ip suppression due to color screening (the lower limit on the 
J ftp formation time) and r = 1.0 fm, a rather conservative estimate of the time needed for the system to reach local 
thermal equilibrium. 

We estimate the maximum energy density in b = U+U collisions with TT orientations. In this case, 
l/S±dN c \ i /dy\ y= () s» 43.9 fm~ 2 while l/S±dN c \ 1 /dy\ y= o ps 31.5 fm -2 in orientation-averaged Au+Au collisions [49]. A 
proposed experimental method for selecting TT orientations gives 1 / S±dN c h/ dy\ y= o « 40 fm -2 although an effective 
trigger has yet to be developed [49] . 

In order to convert energy density to temperature and estimate T/T c , we use lattice QCD results for e(T) from 
Ref. [SO] with T c « 185 MeV. We find, that, for r = 0.25 fm, the temperature exceeds 2.7T C in both U+U and 
Au+Au central collisions. However, for r = 1.0 fm, the temperature is below 2T C in Au+Au collisions while it could 
reach 2T C in central U+U collisions averaged over orientations and may exceed 2T C in TT orientations, see Table [I] 
Consequently, selecting TT configurations in U+U collisions would extend the range of temperatures accessible at 
RHIC to test the hypothesis of sequential quarkonium melting. The temperature in U+U TT collisions would thus 
be high enough to determine whether the T(15) dissociates at T m 2T C [51] , 



Table I: The accessible temperature for r = 1 fm in the 5% most central Au+Au and U+U collisions averaged over spatial 
configurations, ideal TT U+U collisions and those which may be selected experimentally, TT (exp), as described in Ref. [49] . 



System 


T [MeV] 


T/T c 


Au+Au (0-5%) 


351±? 


1.90 


U+U (0-5%) 


366±? 


1.98 


U+U (TT, exp) 
U+U (TT, ideal) 


376 
385 


2.03 
2.08 



V. EXPERIMENTAL SELECTION OF TT AND SS CONFIGURATIONS 



There are several strategies which can be used to select TT and SS (or TT- and SS-enriched) event samples. In 
general, the number of binary collisions in central and mid-central collisions in TT configurations is higher than 
in SS configurations at the same b, leading to higher multiplicities in TT configurations. Therefore TT-enriched 
head-on (b = 0) events could be selected by experimental cuts on the measured charged-track multiplicities. This 
method could be extended to other centralities. The centrality class can be established by employing Zero Degree 
Calorimeters (ZDCs) which measure the energy of spectator nucleons traveling in the forward and backward directions. 
TT-enriched events would be then selected from the centrality class by a multiplicity cut. Results obtained in Monte 
Carlo Glauber simulations show that the charged track multiplicity dN^/drj in central TT collisions is ~ 17% higher 
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than central collisions in the SS configuration at b — 0. The difference decreases with b although diVT r /dr? is still 
10% higher at b = 4.5 fm [46]. Therefore, the multiplicity difference is suitable for offline event selection. Selection of 
Tip and Side orientations in p\J collisions would be probably more challenging because the overall multiplicity will be 
lower. However, the number of binary collisions in p+Tip interactions should be significantly higher than in p+Sidc 
collisions. The relative difference could be large enough for separation to be feasible in those cases. 

Another approach for selecting collision orientation was proposed in Ref. [49]. TT and SS configurations have 
different elliptical eccentricities and, consequently, different magnitudes of anisotropic flow. The 'reduced flow vector' 
was proposed in addition to multiplicity cuts to select SS and TT samples with high purity. 

In any case, selecting a particular collision geometry requires detailed modeling and simulations of experimental 
observables. 

Note that, for CNM studies, high-purity samples are not required: the min-bias data are dominated by SS con- 
figurations. Thus TT-enriched samples with sufficiently longer path lengths L than min-bias collisions in the same 
centrality class would be appropriate for studying nuclear absorption. Experimental cuts could thus be adjusted to 
optimize the balance between the desirable (L) for physics and the efficiency of event selection. 



VI. i?AuAu= A NEW OBSERVABLE FOR J/tp STUDIES 

As shown in Sec. |III| the strength of gluon shadowing is one of the largest uncertainties in the determination of 
cold nuclear matter effects on J/tp production. The effects of shadowing are almost identical for Au+Au and U+U 
collisions averaged over impact parameter. Thus the relative uncertainty can be reduced by studying J/tp production 
in U+U and Au+Au collisions within the same centrality class. As outlined in Sec.|lJ the J ' ftp yield, N^/jf , produced by 
statistical recombination is proportional to the number of charm quarks in the system: N^/J? oc N%, and N C c oc A'bm, 
the number of binary nucleon-nucleon collisions. We define the relative nuclear modification factor -R^uAu as 

ffUU _ dN^/dy / jyAgAu y 
AuA " dN^/dy { AC ) ' { ' 

the ratio of J/tp yields scaled by the number of binary NN collisions, N^ in , in Au+Au and U+U collisions respectively. 

If the J/tp's are produced by coalescence and suppression is indepedent of nuclear absorption and energy density, then 
c-uu ^ i 

n Au+Au ~ i - 

The min bias ratio -RauAu nas an advantage over Rcp, the ratio of the J/tp yield in more central bins relative 
to that in the most peripheral bin, because -RauAu does n0 * depend on shadowing, see Figs. § and [7| As we have 
seen, shadowing can vary significantly between central and peripheral collisions but has a rather weak dependence on 
impact parameter. Indeed, if no selection is made on orientation in U+U collisions, the shadowing dependence will 
cancel in the ratio. 

With the large data samples taken by STAR and PHENIX in 2009 and 2010 (and even larger data sets expected 
for RHIC II), the precision of the J/ip yields will be driven by systematic rather than statistical uncertainties. The 
dominant systematic error (±10%) in pp interactions is the estimate of the integrated luminosity |52) which does not 
cancel between d+Au and Au+Au measurements, as seen in Ref. [36] where the systematic error on i?dAu is much 
larger than that on Rcp- In the case of -RauAuj systematic uncertainties would mostly cancel when Au+Au and U+U 
data are taken by the same detector. 

These features of -RauAu ma ke it useful for testing models of J/tp production and in-medium interactions as described 
in Sec.[T] The dependence of nuclear absorption on b is rather small, but can be detected if large data sets are collected. 
Furthermore R^uAuiv) wm " provide an additional test of models which reproduce RauAu{v)- 



VII. SUMMARY 



We have investigated cold nuclear matter effects on charmonium production in U+U collisions. Such collisions 
provide an interesting opportunity to study J/tp in-medium interactions since model-dependent uncertainties can be 
significantly reduced. We propose a new observable, -RauAu> which is free from some of the uncertainties associated 
with i?AuAu and i?cp- The experimental techniques for quarkonium measurements at RHIC are well established; there 
are abundant Au+Au data; and U+U collisions are planned at RHIC in the near future. Therefore, our proposed 
observable -RauAu represents an additional, important handle on J/tp in-medium interactions. Moreover, the energy 
density achievable in U+U relative to Au+Au collisions is up to « 20% larger, extending the range of energy densities 
available for testing quarkonium suppression due to color screening. 
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Furthermore, U+U collisions would be also useful for cold nuclear matter studies at lower energies, as part of a 
RHIC Beam Energy Scan program. Nuclear absorption is expected to be larger at lower energies [34] while there 
may be antishadowing of the nuclear gluon distribution rather than shadowing. Therefore cold nuclear matter effects 
could be tested more directly by lower energy U+U collisions. 
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